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We predict that graphene is a unique system where disorder-assisted scattering (supercollisions) 
dominates electron-lattice cooling over a wide range of temperatures, up to room temperature. This 
is so because for momentum-conserving electron-phonon scattering the energy transfer per collision 
is severely constrained due to a small Fermi surface size. The characteristic temperature de- 
pendence and power-law cooling dynamics provide clear experimental signatures of this new cooling 
mechanism. The cooling rate can be changed by orders of magnitude by varying the amount of 
disorder which offers means for a variety of new applications that rely on hot-carrier transport. 



A number of interesting and practically useful phenom- 
ena arise when slow cooling between electronic and lat- 
tice systems results in long-lived hot carriers proliferating 
over large spatial scales 1 . Energy transport and energy 
harvesting mediated by hot carriers is utilized in a variety 
of applications (calorimetry, bolometry, infrared and THz 
detectors, etc.). Thermal decoupling of electrons from 
the crystal lattice in most materials takes place at tem- 
peratures of order a few kelvin^. In contrast, the rates 
for electron-lattice cooling in graphene are predicted to 
be very slow in a much wider temperature range [31 
resulting in new optoelectronic and thermoelectric phe- 
nomena as well as other hot-carrier effects [SHS]. 

The inefficiency of the standard cooling pathways me- 
diated by optical and acoustic phonons|3l [4] stems from 
the material properties of graphene. The large value of 
the optical phonon energy, ujq — 200 meV, quenches the 
optical phonon scattering channel below a few hundred 
kelvin[9j; a small Fermi surface and momentum conser- 
vation severely constrain the phase space for acoustic 
phonon scattering |3l |^ . Given the momentum conserv- 
ing character of these processes, which renders them inef- 
ficient, cooling in graphene can be particularly sensitive 
to the effects of disorder. In this Letter, we argue that 
an unconventional, disorder-assisted pathway dominates 
cooling in a wide range of temperatures, explaining key 
features of cooling dynamics observed in recent pump- 
probe measurements ,7J (see Fig. [l]). 

The cooling dynamics reported in Ref. [7j features fairly 
long timescales. The cooling times grow with decreasing 
temperature, from ^ 10 ps at 300 K to ~ 200 ps below 
50 K. This is very different from the dependence ex- 
pected for momentum-conserving scattering by acoustic 
phonons, where the cooling times are predicted to in- 
crease with temperature, reaching a nanosecond scale at 
room temperature for comparable densities [31 3]. The 
observed temperature dependence is also clearly distinct 
from the very steep dependence expected for optical 
phonons, r ~ ghuio/ksT ^ show below, the disorder- 

assisted cooling mechanism yields slow time scales and 
a temperature dependence that closely match the obser- 
vations. In addition, as illustrated in Fig[TJ this mecha- 



(b) 

Decreasing 
T5UB= 300, 200, ISO, 
100, 50, 18 K 




100 200 

Time (ps) 



100 200 300 

Delay Time (pfi) 



FIG. 1: (a) Temperature dynamics obtained from Eq.([T| for 
the lattice temperature values matching those in panel (b). 
Parameter values used: doping n = 50meV and disorder 
mean free path kpl = 20. (b) Carrier dynamics measured 
using the pump-probe technique for varying substrate tem- 
peratures (reproduced from Fig. 2(b) of Ref. 7 ). 

nism explains subtle features such as the prolonged non- 
exponential regime of cooling dynamics and the satura- 
tion of cooling times at low T manifest in the similarity 



between the 50 K and 18 K curves [see Eqs.(|13|),(|l4|]. 

The high impact of disorder on cooling can be under- 
stood by noting that the momentum-conserving acoustic 
phonon processes can only dissipate energy in bits much 
smaller than fc^T. Indeed, since for such processes the 
phonon momenta are limited hy 2kp, the maximal energy 
transfer cannot exceed 2A:bTbg — "ifiskp per scattering 
event (here s and kp are the sound velocity and Fermi 
momentum). The Tbg values are a few kelvin for typ- 
ical carrier densities, i.e. a small fraction of ksT. In 
contrast, disorder-assisted scattering allows for arbitrar- 
ily large phonon recoil momentum values. In this case, 
the entire thermal distribution of phonons can contribute 
to scattering, resulting in the energy dissipated per scat- 
tering of order ksT (supercollisions). This provides a 
dramatic boost to the cooling power. 

For this cooling mechanism, modeling disorder by 
short-range scatterers, we obtain the energy-loss power 



A = 9.62- 



(1) 



where Td (Tph) the electron (lattice) temperature, is 
the density of states at the Fermi level per one spin/valley 
flavor, g is the electron-phonon coupling, and A:BToi(ph) ^ 
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/X. The enhancement factor for the disorder-assisted cool- 
ing over the momentum conserving pathways depends on 
both disorder and temperature: 
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see also Eqs.(12|,(13)]. At room temperature, T^^ph) 



k,i q 

where — vpa ■ k is the massless Dirac Hamiltonian, 
identical for i = 1...7V spin/valley flavors. The electron- 
phonon interaction arises from the deformation potential, 



Hei-ph = X! 5V^(^q + ^-q)"q' 5 = £'/\/2ps2, (4) 
q 

where I? is the deformation potential constant and p is 
the mass density of the graphene sheet. 




300 K, and takmg /x = 100 meV (n ~ 10^^ cm"^) we find 
Tci(ph)/'7BG ~ 50. For kp^ = 20, the enhancement factor 
J/Jo can be as large as 100 times. 

Given the dominance of the disorder- assisted pro- 
cesses, we predict that cooling in graphene is uniquely 
sensitive to disorder. fceTe^ph) ^ A* This sensitivity can 
account for the wide spread of experimentally measured 
cooling times [21 O fT0Hl2] . Slow cooling times arise be- 
cause J scales linearly with the disorder concentration, 
via \/kpt, and with carrier density, n ^ v'^^^l). The 
inverse scaling with kp^ is consistent with the trend of 
cooling becoming faster at higher levels of disorder noted 
in Ref.[TT]. The sensitivity to disorder can be used as a 
knob to engineer cooling rates desirable for specific ap- 
plications. 

The enhancement to phase space may also arise due 
to processes of other types^^. Recently, Castro et. al. 
[m [15] predicted that scattering by flexural phonons can 
dominate the resistivity (momentum relaxation) in free 
standing graphene. In contrast, here we are concerned 
with cooling which is sensitive to both the scattering rate 
as well as the exchange in energy. For flexural phonons, 
we flnd an energy-loss power T dependence resembling 
that in Eq.Q but with a greatly diminished prefactor. 

In our discussion of cooling we shall implicitly assume 
that an effective electronic temperature is established 
quickly via fast carrier-carrier scattering. This is well 
justifled as carrier-carrier scattering occurs on timescales 
of tens of femtoseconds [H [12], far shorter than the 
timescales ^ 1 ps we are concerned with. Below, we an- 
alyze cooling from phonons in the graphene lattice only. 
We note that other phonons (particularly, substrate sur- 
face phonons [T^^ may also contribute to coolingfTT]. 

We proceed to analyze the disorder-assisted cooling 
regime, wherein impurity scattering mediates the ex- 
change of momentum and energy between electron and 
phonon systems. The effect of disorder on electron- 
phonon scattering can be described by the Hamiltonian 

H = ^4,^^o(k)Vk,»+5I^q&q^+^oi-ph+Hdis, (3) where G(p 
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FIG. 2: (a) Kinematics of supercoUisions and normal col- 
lisions at T > Tbg- Phonon momenta are constrained by 
the Fermi surface for normal collisions (white arrows), and 
totally unconstrained for supercoUisions (gph), with the re- 
coil momentum (grccoii) transferred to the lattice via disorder 
scattering or carried away by second phonon. The energy dis- 
sipated in supercoUisions is much greater than that dissipated 
in normal collisions, (b) Feynman diagrams for disorder as- 
sisted electron-phonon scattering processes, corresponding to 
the three terms in Eq.([7|). 

The transition probability for the emission and absorp- 
tion of phonons can be described by Fermi's Golden Rule, 



\^N^5+ + \M4\N^+1)5J^, (5) 



where 5± — i5(ek' ~ T '^q): Q is phonon momentum 
and Ni^^ = 1/(6^"" — 1) is the Bose distribution. In the 
absence of disorder, Eq.(|4| yields the matrix elements 

r(0) 



= .g-y'^'Jk'-kiFq' where the delta function enforces 
momentum conservation. In the presence of disorder, 
possible phonon momenta are unconstrained, taking on 
any value |q| qr ~ ksT/s [see Figj2](a)]. 

We model the disorder potential as a sum of randomly 
positioned impurity potentials, 



Hd,s-5]V'l(r)C/(r)V',(r), {/(r) - ^ F(r 



(6) 



At low disorder concentration, we can describe disorder- 
assisted phonon scattering by dressing the electron- 
phonon vertex with multiple scattering on a single im- 
purity. This gives an expression for the transition matrix 
elements M± which is exact in the impurity potential: 



M± = (k' I GT + tGM^> + fGM^> GT I k) , (7) 



(0) 



— J , , is the electron Green's function, 

T is the T-matrix (scattering operator) for a single impu- 
rity. The three terms in Eq. ([t]) account for the cases when 
impurity scattering occurs before or/and after phonon 
emission (see Fig. 1 (b)). 

As we shall see, the main contribution to cooling will 
arise from phonons with momenta of order qx- Thus we 
anticipate that the virtual electron states, described by 
the Green's functions G(p), are characterized by large 
momenta |p| qx which are much greater than k, k'. In 
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this case, for the off-mass-shell virtual states such that 
"^fIpI ^ k^T, we can approximate G(p) ~ ~TT^^{p)- 
The stiffness of electron dispersion, vp ^ s, along with 
the estimate |p| qt, makes it an accurate approxi- 
mation for all virtual states not too close to the Fermi 
surface. In this limit, as we now show, drastic sim- 
plifications occur because of the particle-hole symmetry 
iJo(-p) = -Hoip). 

We focus on the case of short-range disorder, modeled 
by a delta function potential V{r — r^) = uS{r — rj)(l ± 
az)/2, where the plus (minus) sign correspond to impu- 
rity positions on the A (B) sites of the carbon lattice. 
In this case, a nonzero result for the transition matrix 
element M± is obtained at first order in u. We approxi- 
mate fp,,p = iu(i ± cr^)e*(p'-P)''i + 0[v?) and evaluate 
the first two terms in Eq.([7|. This gives the commutator 
of H^^{c\) and icr^, arising because the virtual electron 
states in the first and second term have momenta p ~ — q 
and p w +q (see above). We obtain 

M± = ^Pf^{\^\{a X q),|k)e»(P'-P±'^)'■^ (8) 

with the phase factor describing the dependence on the 
impurity position. We evaluate the energy-loss power as 

J = > ; I^k',k(ek - ek')/(ek)[l - /(^k')] (9) 



E 

k,k' 



where /(e) = l/(e'^(<=-'') + 1) are Fermi functions, Wk',k 
is the transition probability, and ek — ek' is the energy 
exchanged in each scattering event. 

In the degenerate limit, ksT <^ ii, the sum over k and 
k' is conveniently factored into separate integration over 
energies and angles X^k.k' ^ ['^(m)]^ //c^ede' // 
One of the energy integrals is eliminated by the delta 
functions 5(ek' — ek i'^q)- The second energy integral is 
evaluated using the identity de/(e)(l — f{e + ujq)) ~ 

^1 [-^iSq + l] ' where is the Bose distribution function 
evaluated at the electron temperature. With the electron- 
phonon matrix element given by Eq.(|8]), and using the 



angle-averaged quantity (|(k'|((T x q)z|k) 
we obtain an expression 
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where no is impurity concentration. Integration yields 
Eq.([T]), where we used an expression for the mean free 
path kpt = 2fL^vl./{u'^no)^. 

We can make a comparison with the normal (momen- 
tum conserving) processes[31l3], where the cooling power 
is Jo — B{Tei — Tph), where B = nN\hi'{fi)kps'^kB, and 
A = g^v^fi) is the dimensionless electron-phonon cou- 
pling. Linearizing Eq.([T]) in AT = Td — Tph, we find that 
this contribution dominates over Jo ^-t temperatures 
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Taking kpl = 20 for a rough estimate, we see that the dis- 
order assisted cooling channel dominates for T ^ STbg- 
The crossover temperature can be controlled by gate volt- 
age, since Tbg oc y/n. For typical carrier densities n this 
gives a crossover temperature T* of a few tens of kelvin. 

Interestingly, the cooling times describing relaxation 
to equihbrium, ATei(t) = e"(*"*'')/^AToi,o, exhibit a 
nonmonotonic T dependence for T ^ T^,. Accounting 
for both the disorder-assisted and momentum-conserving 
cooling, the relaxation dynamics can be described as 



dQ/dt = - J - Jo, 



(12) 



where Q is the electron energy density. Taking Q = 
CAT, with C = aT the heat capacity of the degenerate 
electronic system (a — ^Nv{ij)k'^), we find 



= — T 



(13) 



B 

T a aT 
The cooling time increases with T at T < T* and de- 
creases at T > T*, reaching maximal value at T = T*. 
The non-monotonic temperature dependence provides a 
clear experimental signature of the competition between 
different cooling pathways. 

To describe the cooling dynamics both near and away 
from equilibrium, we used non-linearized quantities, Q — 
laT^^j and Eq.(|l|), with the deformation potential con- 
stant D — 20eV|31[l], the electron temperature initial 
value Toi^o = 3 • 10^ K ~ wq and other parameter val- 
ues cited in Fig{T] caption. For the parameters used, 
~ 15 K. The resulting dynamics, shown in Figjl|a), 
reproduces the main features seen in the data. We note 
that the long time behavior is insensitive to the choice of 
Toi^o; only the dynamics at short times are affected. 

The non-exponential behavior seen in the data at short 
times can be understood by analyzing the regime T^i S> 
Tph. Approximating J « AT^^ and suppressing Jp, we 
obtain a l/(t — to) dynamics: 

7ci,o 



T,i{t) 



l + iA/a)it-to)T,i^o' 



(14) 



The dynamics at intermediate times, where Tci ?i ?ph. 



can be found by directly solving Eq.(12). We obtain 
--{t^tn) = F{T,i{t)/Tph) - T(Tei,o/rph), (15) 

T 

where F(x) = 2^3 arctan[(l + 2a;)/^/3] -ln[(a;3 - l)/(a;- 
1)'^] (we suppressed the Jo term which becomes impor- 
tant only for T T^, and only at long times). This so- 
lution, with T taken from Eq.(13), was used to generate 
Fig. [l](a), yielding results strikingly similar to the data. 

Interestingly, both at T < T* and T > T* the super- 
collision frequency remains lower than that for normal 
processes, (Wq) = ^XknT^]^. We can define the av- 
erage collision frequency as 

Ek ,k'W^k^k/(£k)[l-/(ekO] 
Ek/(ek)[l-/(ek)] 



{W) = 



(16) 
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Evaluating the integrals as above and setting Tci = ?ph, 
we find {W) = j^gXkBTln^ < (Wq). The low value 
for (W) means that normal collisions produce the domi- 
nant contribution to resistivity even when their contribu- 
tion to cooling is totally overwhelmed by supercollisions. 

Finally, we analyze cooling in free standing graphene 
in the absence of disorder. In this case, an important 
contribution arises due to flexural phonons |14lfT5] . which 
contribute to the deformation tensor via Uij = l/2(diUj + 
djUi+dihdjh), with u and h the in-plane and out-of-plane 
displacements. Flexural modes have quadratic dispersion 
(Iiq = K|q|2 with Ksa 4.6 X 10""^ m^s-^tll [H]. Electron- 
phonon coupling is described by the same deformation 
potential as above, Eq.Q. 

The processes involving pairs of nearly counterprop- 
agating flexural phonons are analyzed as follows. Us- 
ing the momentum representation, hq — ^Jfij2pCjq{hq + 

&lq), we consider the emission/absorption of two flexu- 
ral phonons with momenta qi and q2. For T ^ T^'q = 
hnk^p (for typical densities, Tgg' is well below IK), we 
can set qi ~ — q2 = q, yielding the transition probability 
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k',k- 



(17) 

5± = (5(ek — Ck' i 2wq). Here the matrix element is M — 
^ [aiia and the coherence factor is |(k'|k)|2 = [l ± 

cos(0k — ^k')] /2j with the plus (minus) sign for intra-band 
(inter-band) processes. This gives the energy-loss power 



9 



i2hu,y {N^+ l)^iV°' - < {N^l+ 1) 



where - = Ht^)] ' I^- We note that the 

above expression vanishes when Td = Tph- 

We linearize jTflox in AT = T^i — Tph to obtain 



Jfle. = AiT^AT, Ai = 0.12 



ND^v^{ti)kl 



(18) 



which scales with T the same way as Eq.Q. Flexural 
phonons dominate over the one-phonon contribution at 
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(19) 



The value T^°^ is considerably larger than for 
disorder-assisted cooling estimated above. A comparison 
with EqQ yields Jncx/J ~ kpl/^OO which is small for 
typical kpi.Thus the contribution ([l8| is relatively weak 
under realistic conditions. For graphene on substrate this 
contribution is further diminished as flexural modes get 
pinned, gapped, and stiffened by the substrate. 

Besides carrier dynamics, cooling can also be probed 
by transport measurements through bolometry |19j . de- 
scribed by the thermal impedance, i?th = {dP/dAT)~^. 



Here V is the power pumped into the system (say, via 
Joule heating), which is balanced by J" in a steady state. 
Hence, i?th temperature dependence can be used as a 
diagnostic for the processes dominant in cooling. For 
disorder-assisted cooling, extending Eq.([T]) to ^ < fceT, 
we approximate R^^ = AT^(l -I- c(fcBT/^)^) for mono- 
layer graphene, with c a constant of order unity. This 
temperature dependence is markedly different from i?th 
predicted for momentum conserving channels [20] . 

In summary, graphene stands out as a unique system 
where disorder-limited cooling is the leading contribu- 
tion over a wide range of temperatures, including room 
temperature. As a result, varying the amount of disor- 
der can be used as a knob to tailor and control a va- 
riety of hot carrier effects in graphene. Tuning disor- 
der can be achieved by well established techniques, in- 
cluding current-annealing and using different substrates 
(e.g. Si02 or BN). The characteristic dependence, 
Eq.([T]), and power-law coohng dynamics, Eq.( [T4|, m akes 
this new regime easy to identify in experiments |21l I22j . 
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